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The quantum phase transition from a spin-Peierls phase with a small Fermi surface to a para- 
magnetic Luttinger-liquid phase with a large Fermi surface is studied in the framework of a one- 
dimensional Kondo-Heisenberg model that consists of an electron gas away from half filling, coupled 
to a spin- 1/2 chain by Kondo interactions. The Kondo spins are further coupled to each other 
with isotropic nearest-neighbor and next-nearest-neighbor antiferromagnetic Heisenberg interactions 
which are tuned to the Majumdar-Ghosh point. Focusing on three-eighths filling and using the 
density-matrix renormalization-group (DMRG) method, we show that the zero-temperature tran- 
sition between the phases with small and large Fermi momenta appears continuous, and involves 
a new intermediate phase where the Fermi surface is not well defined. The intermediate phase 
is spin gapped and has Kondo-spin correlations that show incommensurate modulations. Our re- 
sults appear incompatible with the local picture for the quantum phase transition in heavy fermion 
compounds, which predicts an abrupt change in the size of the Fermi momentum. 

PACS numbers: 75. 30. Mb, 71.10.Hf, 71.10.pm, 71.27.-|-a 
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I. INTRODUCTION 

The Kondo lattice model (KLM) was historically intro- 
duced to describe the competition between singlet forma- 
tion and magnetic ordering in heavy fermion systems.^ In 
heavy fermion materials, localized /-shell electrons hy- 
bridize with itinerant electrons. Depending on whether 
the / electrons participate in the formation of the Fermi 
surface (FS) or not, the latter may be large or smalL^ 
The Fermi momentum kp is large in the sense that the 
FS encloses a volume that counts both the number of con- 
duction electrons and local moments. ^''^ This is in con- 
trast to a small Fermi momentum kp, whose FS encloses 
a volume that counts the number of conduction electrons 
only. 

An earlier issue debated on the KLM was whether the 
formation of a large FS is consistent with Luttinger's 
theorem.^ In other terms, the question was whether the 
KLM can account for a large FS given that the / elec- 
trons are represented by their spin degrees of freedom 
only. Investigations of the one-dimensional (ID) KLM 
through numerical-^" and analytical^''^- approaches have 
yielded rather consistent evidence for a ground state with 
a large FS in the paramagnetic regions of the model. 
Yet recent numerical calculations have suggested the ex- 
istence of a second phase where the FS is small. Such 
a phase will necessary have a broken-symmetry ground 
state, otherwise it would be inconsistent with Luttinger's 
theorem.liS For reviews on the ID KLM, we refer the 
reader to Refs. d andlT^ 

In the last decade, the question of the size of the 
FS in the KLM has gained renewed interest in connec- 
tion with quantum criticality and the related non-Fermi- 
liquid phases of heavy fermion materials. The local pic- 
ture for the quantum phase transition (QPT) in these 
compounds predicts that the size of the FS would change 
abruptly at the quantum critical point J^i^ The com- 



posite quasiparticles forming the large FS are projected 
to breakdown as the system is driven across the critical 
point, leaving behind a small Fermi volume that counts 
the number of conduction electrons only. It is still un- 
clear how such a sudden change of the FS is consistent 
with a second-order phase transition. 

One-dimensional models, for which powerful methods 
of solution are available, are currently the primary tool 
for gaining reliable information about QPT in the KLM. 
However, magnetic orderings that break spin-rotational 
symmetry are prohibited in ID. Hence, it is necessary 
to study the QPT between the paramagnetic phase with 
a large FS and alternative phases of the Ising or spin- 
Peierls type.— If the / electrons are in an Ising or a 
dimerized phase, they would remain decoupled from the 
conduction electrons at low energies also in the presence 
of a sufficiently small Kondo coupling as compared to the 
gap. Consequently, the FS would be small. A QPT to- 
ward a ground state with a large FS would occur upon 
increasing the Kondo coupling. This latter phase is pre- 
sumably a Luttinger liquid (LL) or a spin-gapped phase. 

In this paper, we study the evolution of the ground 
state of a ID KLM from a spin-Peierls phase with a small 
Fermi momentum kp to a LL phase with a large Fermi 
momentum kp. To this end, we augment the conven- 
tional KLM with isotropic nearest-neighbor and next- 
nearest-neighbor spin-exchange interactions among the 
Kondo spins, and tune them to the Majumdar-Ghosh 
point The inclusion of next-nearest-neighbor coupling 
is a crucial ingredient of our study, as it supports a 
broken-symmetry ground state with a small FS for small 
Kondo couplings. 

Focusing on three-eighths filling and using the density- 
matrix renormalization group method (DMRG)fii we 
find a zero-temperature transition that is more com- 
plex than the predictions of the local critical theory. 
In particular, we identify an intermediate spin-gapped 
phase in between the spin-Peierls and LL phases where 



2 



the Fermi momentum cannot be defined. Instead, the 
electron momentum distribution function n(fc) displays 
a shallow peak at a new characteristic momentum k* 
that lies in between fcf. and kp, and which shifts to- 
ward kp on going from the spin-Peierls to the LL phase. 
Concomitantly, there is a maximum at 2k* in the mag- 
netic structure factor S{k) of the Kondo spins. We 
show that this behavior can be understood as a conse- 
quence of the magnetic frustration induced by compet- 
ing spin-exchange couplings generated by the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction. In contrast 
to n{k) and S{k), neither the Fourier transform of the 
local conduction-electron density nr{k) nor the density- 
density correlation function C{k) show any special sig- 
natures related to k* . Rather, the transition from the 
spin-Peierls to the spin-gapped phase is manifest in C{k) 
by the smearing of a cusp at 2kp and the emergence of a 
peak at 2kp. Our results appear incompatible with the 
local picture for the QPT in heavy fermion compounds. 

The remainder of the paper is organized as follows. 
In Sec. |TT] we present the KLM under study, along with 
details of our DMRG code. A comprehensive set of re- 
sults for the dimer order parameter, the spin velocity, the 
electron momentum-distribution function, the magnetic 
structure factor of the Kondo spins, and various density 
correlations are described and analyzed in Sec. IIIII We 
conclude in Sec. IIVI with a discussion of the resulting 
phase diagram and its relevance to heavy fermion com- 
pounds. 

II. THE MODEL 

In this paper we study the following KLM 

L-l L 

H = -tJ2 {4.aC^+l.a + H.C.} + Jk ^ S, ■ n 

i—l,a i—1 
L-1 L-2 
+ Jhi ^ Si ■ Si+i + Jh2 ^ Si ■ Si+2, (1) 
1=1 i=l 

describing a ID tight-binding conduction band with 
the hopping term t, interacting via an on-site spin- 
exchange (Kondo) interaction Jk with an array of lo- 
calized spins. Here, ^ creates a conduction electron 
with spin-projection a at site i. Si represents the local- 
ized spin at site i, and Ti = ■^J^a a' cl o-Ci.o-' (t?)^,^' is the 
conduction-electron spin density at that site. The con- 
duction electrons and spins reside on an L-site lattice 
with open boundary conditions (OBC). In addition to 
the Kondo interaction, the localized spins interact among 
themselves via the nearest-neighbor and next-nearest- 
neighbor Heisenberg spin-exchange terms Jhi and Jh2, 
respectively. To avoid the onset of ferromagnetism^ we 
set Jhi = t/2, while Jh2 = J Hi/ '2 is tuned to the well- 
known Majumdar-Ghosh point,— whose corresponding 
ground state of the decoupled spin chain is a perfect 
dimerized state (for even L). 



As emphasized above, the inclusion of Jh2 = Jhi/'^ > 
is a crucial difference from previous DMRG studies of 
the KLM^^^^^^^^^ This additional frustration opens a 
gap in the spectrum of the isolated spin chaiuf^i enabling 
the study of the transition from the broken-symmetry 
dimerized phase with a small FS to the LL phase with 
a large FS. From a technical standpoint, the dimeriza- 
tion gap significantly reduces the numerical effort that 
is needed to obtain reliable results as compared to the 
case where Jh^ = Jh2 = 0,—"— due to the short-range 
spin-spin correlations that develop. 

We computed the ground state of the Hamiltonian of 
Eq. (P) using the DMRG method with OBC. We retained 
between 256 and 512 states in the two external blocks, 
keeping track of the total number of electrons and the z 
component of the total spin projection Sz as good quan- 
tum numbers. The maximal truncation error was in the 
order of lO"** when 256 states where kept and in the 
order of 10~^ when 512 states where kept. We studied 
different lattice sizes up to L = 64 sites with Jk varied 
in the range < JK/t < 16. For concreteness we set 
the conduction-electron filling equal to n = 0.75, which 
is close to but off half filling, and is rather convenient to 
tackle numerically. We briefly comment on other filling 
factors at the end of the paper. All results presented be- 
low are restricted to zero temperature. The lattice size 
is L = 64 unless stated otherwise. 



III. RESULTS 

It is instructive to consider first the limits of small and 
large Jx, where different phases are expected. When 
Jif = 0, the conduction electrons and spins are decou- 
pled, forming independent chains. The spin chain, be- 
ing tuned to the dimerized Majumdar-Ghosh phase, is 
gapped due to the breaking of translational symmetry. 
The electron chain is gapless in both the spin and charge 
sectors, as is the overall system. Due to the gap in 
the spectrum of the isolated spin chain, the Majumdar- 
Ghosh phase is expected to be stable against the inclusion 
of a small Jk- 

In the opposite limit Jk ^ Jhi , the conduction elec- 
trons and spins bind to form localized singlets. At tem- 
peratures below Jk there is no thermal energy to break 
the Kondo singlets, which can be viewed as holes in the 
underlying Kondo spin texture. Hence, by analogy with 
the case where Jh2 = the system is described by the 
t ~ Ji — J2 model for holes, with t — > t/2, Ji = Jhi, 
and J2 — Jh2- Note that the original electronic filling 
factor n (assumed to be smaller than 1) is converted by 
this mapping to the hole filling factor nhoie = The 
corresponding ground state of the t~ Ji~ J2 model is ex- 
pected to be a paramagnetic LL for rihoic = 0.75, which 
differs in symmetry from the dimerized phase at small 
Jk- Thus, a QPT should occur upon increasing Jk- Be- 
low we confirm this scenario and thoroughly discuss the 
nature of the phase transition. 
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FIG. 1: (Color online) (a) The dimer order parameter \D\ of 
Eq. plotted as a function of J/c/t for n — 0.75 and zero 
temperature. For Jk = 0, the system is in a perfect dimer 
state, corresponding to \D\ — 0.25. As the Kondo interaction 
Jk is switched on the dimerization progressively decays until 
it vanishes for Jx/i ^ 1-7. Note the particularly sharp slope 
around JK/t = 1.3-1.4. (b) The spin velocity Vs of Eq. ^ as 
a function of Jk /t, plotted for different lattice sizes L and n — 
0.75. Different qualitative behaviors are observed for weak, 
intermediate, and strong couplings. In particular, a spin gap 
opens in the intermediate-coupling regime 1.3 < JK/t ^ 4. 

A. Dimer order parameter 

The first quantity we study is the dimerization order 
parameter D, defined as 

D^^J2(^^■s.+l){~ly■ (2) 

i—l 

Here the alternating (—1)* factor comes to distinguish 
the dimerized phase from a translational-invariant state. 
Figure [Ija) depicts the evolution of |D| with increasing 
Jk- For Jk — there is perfect dimerization, corre- 
sponding to \D\ = 0.25. This value of |D| stems from 
the fact that each spin Si forms a perfect singlet with 
one of its neighbors and is uncorrelated with its other 
neighbor. Consequently, {Si ■ Si+i) equals —3/4 (0) for 
odd (even) i. With increasing Jk > 0, the dimeriza- 
tion order parameter decreases first gradually and then 
sharply around JK/t = 1.3—1.4. The sharp slope in the 
latter regime suggests a rather rapid change in the na- 
ture of the ground state. Eventually \D\ vanishes above 
JK/t w 1.7, indicating the loss of any remnant of the 
dimer state that forms at small Jk- It should be noted 
that D shows no significant size dependence due to the 
short-range spin-spin correlations that are involved. 

B. Spin velocity 

When JK/t 1-7 translational symmetry is restored, 
as signaled by the vanishing of 13. In order to better 
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FIG. 2: (Color online) The spin gap As{L) as a function of 
1/L, plotted for n = 0.75 and J/c/i = 0.5, 2.5, and 4. The 
solid lines are linear extrapolations to L — > cxd. In contrast to 
JK/t — 0.5 and 4, where As{L) extrapolates to zero, a finite 
spin gap As{L -> oo) ~ 0.077t is found for JK/t = 2.5.— 



understand the nature of this new phase we computed 
the spin velocity v^, defined as 

VsiL) ^ As{L)L. (3) 

Here As{L) is the elementary singlet-triplet excitation 
energy for a system of size L. Figure [IJb) shows Vs{L) 
for different lattice sizes as a function of Jk- When 
JK/t ^ 1-3, Vs{L) strongly depends on Jk, reflecting 
the progressive formation of a composite quasiparticle 
made up of the local spins and the itinerant electrons. 
By contrast, Vs{L) is almost independent of both Jk and 
L when 4 < Jk /i- This behavior can be understood from 
the fact that the system is rather well described in this 
regime by the < — Ji — J2 model discussed above, which 
forms a LL when rthoio = 0.75. In the intermediate range, 
1-3 ^ >/k/^ ^ 4, Vs{L) depends strongly on both Jk and 
L, indicating the emergence of a nonzero spin gap As for 
L oo- 

To support this interpretation we have plotted As(L) 
vs. L in Fig. [2l for three representative values of JK/t- 
For both small and large Jk (represented by Ja'/^ = 0.5 
and 4, respectively), As{L) extrapolates nicely to zero as 
L — > cxD, indicative of a gapless state. However, for the 
intermediate value of JK/t = 2.5, As(L) extrapolates 
to the finite spin gap Ag/t « 0.077 as L — >■ oo.2S Such a 
global spin gap is neither consistent with a spin-dimerized 
phase nor with a LL phase. A similar spin gap was found 
throughout the range 1-3 < /t < 4, though the precise 
boundaries of this new phase are somewhat difficult to 
pin downi^ 



C. Electron momentum-distribution function 

Next we address the size of the FS and its evolution 
upon going from Jk = to JK/t = 16, thereby cross- 
ing the three different regimes of small, intermediate and 
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FIG. 3: (Color online) The electron momentum-distribution 
function of Eq. Q, for n — 0.75 and different Kondo couplings 
JK/t = 0,1, 1.3, 2, 3, 4 and 16. When Jk = 0, the con- 
duction electrons form a free band, whose exact momentum- 
distribution function is depicted by the red dots. Accord- 
ingly, there is a sharp step in the momentum-distribution 
function at the (small) Fermi momentum kp — Tm/2 — 
0.3757r. The step at kp persists for small Kondo interactions, 
Jif/t ^ 1.3. In the opposite limit of large Jk (represented by 
JK/t = 16), there is a new step at the large Fermi momentum 
kp = 7r(n + l)/2 = 0.8757r. The corresponding momentum- 
distribution function compares rather well with that of a non- 
interacting tight-binding chain of equal length and the filling 
n' = 1 + n — 1.75, scaled down by a factor of n/n' — 3/7 
(black dashed line). In contrast to the limits of small and 
large Jk, there is no clear sign of a FS for the intermediate 
couplings JK/t = 2 and 3. 



large Jk signaled by D and Vs in Figs. [T] and [D As 
pointed out earlier in the introduction, previous com- 
putations^ on the KLM with Jh2 = concluded that a 
large FS forms in the absence of symmetry breaking. Due 
to the gap in the spectrum of the dimerized Majumdar- 
Ghosh phase, the conduction electrons and spins remain 
decoupled at low energies even in the presence of a small 
Jk, thus forming a small FS. On the other hand, a large 
FS should be recovered when Jk ^ t, Jh^ , as the result- 
ing behavior should basically reproduce that of Jh2 = 0- 
In order to study the transition between these vastly 
different Fermi momenta, we computed the electron 
momentum-distribution function, defined as 



(4) 



i=l 



This definition of n(fc) differs from the conventional one 
n{k) = EAc 



Ck,, 



due to the OBC used. Nevertheless, 
it contains similar information on the FS, as we show be- 
low. In particular, the two definitions must coincide in 
the thermodynamic limit L — )■ cx), provided that transla- 
tional symmetry is not broken. To see this we note that 
the correlator (c| ^C/,/2,cr) with fixed i — L/2 becomes 
independent of the boundary conditions for L ^ oo. 



Converting to periodic boundary conditions and in the 
presence of translational invariance, Eq. ^ reduces then 

^iEct {ci,cr(^k,a) + (cik^^c^kM) , which is nothing but 
the conventional momentum-distribution function. Here 
we made use of the fact that (cj^ ^Ck^a) and (c^j, jjC_fc o-) 
are identical by virtue of inversion symmetry. 

Figure [3] shows n{k) for different Kondo interactions 
at the fixed filling factor n — 0.75. For Jk = 0, one can 
compute n{k) exactly from the single-particle eigenstates 
of the decoupled conduction-electron chain. The exact 
results, depicted by the red circles, essentially coincide 
with our DMRG data, serving as a critical check for the 
accuracy of our code. As expected, there is a sharp step 
in n{k) at the small Fermi momentum kf, = 7rn/2 = 
0.3757r of the free band. Note that n{k) oscillates as a 
function of k, and can become both negative and may 
exceed one. These finite-size effects are eliminated in the 
thermodynamic limit L — )■ oo. We emphasize, however, 
that the definition of Eq. does not require that < 
n{k) < 1. 

As soon as the Kondo interaction is switched on, the 
system gradually looses the sharpness of the FS, as ex- 
pected of an interacting ID system. Nevertheless, a 
clear Fermi momentum can still be observed at kp for 
Jx/t ^ 1.3, i.e., in the range where a sizeable dimerized 
order persists (see Fig. [IJ. In the strong-coupling limit 
4 < JK/t (note that 4t is the free conduction-electron 
band- width), a new well-defined Fermi momentum ap- 
pears, this time at kp — 7r(7i + l)/2 = 0.8757r. This value 
of kp corresponds to a FS which encloses a volume that 
counts both the number of conduction electrons (filling 
factor of n = 0.75) and the number of local moments 
("filling factor" of n = 1). For comparison, the black 
dashed line represents the exact momentum-distribution 
function for a noninteracting tight-binding chain of equal 
length and the filling n' = 1 -I- n = 1.75, scaled down by 
a factor of n/n' — 3/7. This latter renormalization re- 
fiects the fact that the actual conduction-electron filling 
in our system [corresponding to the integrated weight of 
n(k)] is n rather than n' ^ 1 + n. The agreement is quite 
surprising. 

In contrast to the weak- and strong-coupling regimes, 
there is no well-defined Fermi momentum for the inter- 
mediate couplings Jif/i = 2 and 3. Instead, the sharp 
structure at kp is rapidly smoothed and suppressed, and 
a new feature appears at a characteristic momentum k* 
located in between kp and kp. This new feature, whose 
position is marked by the arrows in Fig. |31 first ap- 
pears for Jif/t ~ 1.3 as a small and shallow peak at 
k* = 0.6257r. It continuously shifts toward kp upon in- 
creasing Jk, until it coincides with the new Fermi mo- 
mentum kp for Jif/t = 4. As seen from the Friedel 
oscillations and the density-density correlation function 
presented below, k* is not associated with a new Fermi 
momentum. Rather, it reflects the spin-spin correlations 
that develop in this range due to the combination of Jh^ , 
Jh2 , and the RKKY interaction mediated by the conduc- 
tion electrons. In other words, k* stems from the back 
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FIG. 4: (Color online) The magnetic structure factor of 
Eq. ((5}, plotted for n = 0.75 and the same values of Jk 
as in Fig. [3l For Jk = 0, the Kondo spins are locked in a 
perfect dimer state. The exact magnetic structure factor in 
this case (red dots) has a broad peak at fc = tt. Once Jk is 
switched on, a sharp cusp gradually develops at 2kp — 0.757r, 
reaching a maximum for JK/t ~ 1.3. Upon further increas- 
ing Jk, the cusp smoothens and shifts continuously all the 
way to g = 27r — 2kp = 0.257r, where a new cusp forms at 
large Jk- In the limit of large Jk, represented by Jx/i = 16, 
S{k) well matches the structure factor of the corresponding 
t — J\ — J2 model, which is shown for comparison by the black 
dashed line. Inset: Evolution of the cusp at 2k% as it first 
begins to move. Here crosses, pluses, stars, circles, squares, 
and triangles correspond to Jjf/t = 1-1, 1-2, 1.3, 1.4, 1.5, 
and 1.6, respectively. For Jx/i < 1-3 (blue lines), the cusp is 
pinned to 2k%, growing in magnitude with increasing Jk- For 
1.4 < Jx/i (red lines), the cusp is shifted to incommensurate 
modulations. 



action of the presumably frustrated spins on the conduc- 
tion electrons. 



D. Magnetic structure factor 

The structure that n{k) develops at k* in the interme- 
diate phase is rather small. We now show that a much 
more pronounced feature is found at 2k* in the magnetic 
structure factor of the Kondo spins, defined as 



S{k) = Y,{S^ ■ Sl/2) cos[fc(z - i/2)]. 



(5) 



Similar to n{k), Eq. ([5]) is defined to match the OBC 
used. It reduces in the thermodynamic limit and in 
the absence of translational symmetry breaking to the 
Fourier transform of the spin-spin correlation function 

Figure |4] shows the evolution of S{k) with increasing 
Jk- For — 0, when the spin chain locks in a perfect 
dimerized ground state, S{k) equals 0.75[1— cos(fc)]. This 
curve, depicted by the red dots in Fig.|31 has a maximum 
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FIG. 5: (Color online) The magnetic structure factor of the ef- 
fective Heisenberg Hamiltonian with different spin-exchange 
couplings Ji, J2, - - - , J(j, up to a distance of 6 lattice sites. 
Starting from the Majumdar-Ghosh point and gradually in- 
creasing the relative strengths of J2 to Je, we are able to shift 
the main modulation in S{k) from fc = tt to 2fcf' = 0.757r, and 
then all the way to q = 0.257r, in close analogy to the peak 
positions observed in Fig. [3] The incommensurate modula- 
tions that are seen in Fig. |4] can thus be understood as due to 
the frustrated RKKY interaction, whose magnitude increases 
with Jk- 



at fc = TT and is well reproduced by our DMRG data. In 
the opposite limit of a large Jk, S{k) reduces to the 
magnetic structure factor of the t — Ji — J2 model, which 
has a cusp at g = 27r — 2kp — 7r(l — n) = 0.257r [note 
that q and 2kp are equivalent by virtue of S{k) = 5'(27r — 
fc), and will be regarded as synonymous hereafter]. The 
transition between these two limits proceeds as follows. 

When Jk is switched on, the spins couple to the con- 
duction electrons. As a result, S{k) gradually deforms 
from its perfect dimerized profile and develops a sharp 
cusp at 2fc|^ — 0.757r, as can be seen for J^/t = 1 and 
1.3. The evolution of this cusp is tracked in the inset. 
With increasing Jk the cusp initially grows sharper, un- 
til reaching a maximum for JK/t ~ 1.3. Upon further 
increasing Jk (1.4 < JK/t < 4), the cusp smoothens 
and shifts continuously all the way to q — 0.257r. In- 
terestingly, the cusp at 2kp first begins to shift at the 
same value of Jk where the dimerized order parameter 
D acquires its sharpest slope. Furthermore, in the in- 
termediate regime 1.4 < JK/t < 4, the associated peak 
in iS'(fc) occurs at q* — 2tt — 2k* , where fc* is the posi- 
tion of the additional structure found in n{k)- Hence the 
two features are intimately related. Lastly we note that 
the local maximum that is seen at fc = tt for JK/t = 2 
is a finite-size effect that decreases in magnitude with 
increasing system size. This is in contrast to the incom- 
mensurate peaks at q* which depend only weakly on L. 

In order to understand the origin of the incommensu- 
rate modulations in S(k), we studied an effective spin 
Hamiltonian which comes to mimic the RKKY interac- 
tions present in the KLM. Specifically, we consider a frus- 
trated Heisenberg model with different spin-exchange in- 
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teractions Ji , J2 , . . . , Je, up to a distance of 6 lattice sites. 
Such an effective Hamiltonian with adequate couplings is 
expected to properly describe the spin dynamics in the 
KLM up to moderately large values of Jk- This descrip- 
tion must eventually break down for large Jk, when the 
conduction electrons and spins tightly bind to form lo- 
calized singlets. 

Figure [5] shows the magnetic structure factor S{k) for 
the effective Heisenberg model. As can be seen, we are 
able to generate incommensurate modulations in S{k) 
by varying the relative strengths of the different spin- 
exchange interactions. In particular, starting from the 
Majumdar- Ghosh point and gradually increasing the rel- 
ative strengths of J2 to Jg, we are able to shift the main 
modulation from k = tt (perfect dimers) to 2kp — 0.757r, 
and then all the way to q — 0.257r. Note that the peak 
positions in Fig. [5] are quite similar to those in Fig. 21 
though their profiles progressively deviate from those in 
Fig. 21 upon increasing the couplings. This is to be ex- 
pected from the approach to strong coupling, where the 
KLM reduces to the t — Ji — J2 model. These results in- 
dicate that the effective RKKY interaction between the 
Kondo spins is responsible for the incommensurate mod- 
ulations observed in Fig. SI 



E. Density correlations and charge structure factor 

In order to further prove the spin origin of the feature 
that n{k) displays at fc*, we proceed to show results on 
density correlations and the charge structure factor. Fig- 
ure [5] depicts the smoothed Fourier transform of the local 
density nr{k), defined as 



q=27t-2kp 



nr{k) = 



W{i) cos 



1 



— 0.75 is the filling factor, hi equals 



(6) 



Here, 

with hi^a- — cl^Ci^cr, and W{i) is a smooth windowing 
function introduced to avoid spurious edge effects. ^^'^^ 
The additional term I/J with / = J2f=i{'''^i ~ n)W{i) 
and J — J2f=i ^ii) comes to correct for the weighted 
average of fii, which slightly deviates from n. Its role is 
to remove an artificial feature near k — introduced by 
W{i). 

The Fourier transform of the local density oscillations 
is a measure of the Friedel oscillations that develop due 
to the OBC used. It thus offers a direct way to probe the 
Fermi momentum, which is manifest as a peak at 2kp 
and its higher harmonics Figure |6] shows the Friedel 
oscillations for three values of Jx/t — 0, 1.3 and 16. For 
Jk = we observe a single peak at 2kp — O.TSvr, which 
well agrees with an exact evaluation of Eq. © using the 
single-particle eigen-modes of the decoupled conduction- 
electron chain (red circles). In the opposite limit of large 
Jk, there is one peak a.t q = 2tt — 2kp = 0.257r and 
another peak at 2q = O.Stti^ Both structures agree fa- 




FIG. 6: (Color online) The smoothed Fourier transform of 
the local density defined in Eq. for n = 0.75 and diflferent 
Kondo interactions. For Jk ~ 0, there are pronounced Friedel 
oscillations at the modulation wavelength 2kp = 0.75n, in 
good agreement with an exact evaluation of Eq. ® using the 
single-particle eigen-modes of the isolated tight-binding chain 
(red circles). At large Jk, represented by JK/t — 16, there 
are Friedel oscillations at both g = 27r — 2kp = 0.257r and 
2q = O.Stt, in accordance with |nr(fc)|^ for the corresponding 
t~Ji — J2 model (black dashed line). Surprisingly, for Jk /t = 
1.3 there are simultaneous modulations at 2kp, q, and 2q, as 
if there were two coexisting Fermi momenta kp and kp (note 
that 2q — 2tv — 4fcf = 47r — Akp corresponds to both 4fef and 
Akj^). 



vorably with similar calculations for the corresponding 
t — Ji — J2 modelf2^ whose results are shown by the black 
dashed line. Remarkably, for the intermediate coupling 
JK/t — 1.3 we find modulations at all three momenta 
2kp, q, and 2q, as if there were two coexisting Fermi 
momenta^^ kp and kp (note that 2q corresponds both 
to 4:kp and 4fc|^). These results are in stark contrast to 
the electron momentum-distribution function of Fig. [31 
which shows only a single step at kp for this value of Jk- 

As a function of Jk, the peak at 2k p persists from 
Jk = up to JK/t ~ 1.7, which is the point where the 
dimerized order parameter \D\ first vanishes (see Fig.[T]). 
Above JK/t ~ 1.7 there are no discernable signatures 
left at 2k p. The peak at q appears continuously at small 
Jk, and persists all the way up to large Jk- It evolves, 
however, in a nonmonotonous fashion, first growing in 
magnitude before decreasing toward its asymptotic t — 
Ji — J2-niodel form. Interestingly, no special signature 
appears at the incommensurate momentum q* = 27r — 
2k*, indicating that k* is unrelated to the FS. 

Finally, in Fig. [7] we show the Fourier transform of the 
conduction-electron density-density correlation function, 
defined as 

L 

C{k) = Y,{{n, - n){nLi2 - n)) cos[k{i - L/2)]. (7) 

i=l 

Similar to previous plots, the DMRG data well repro- 
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FIG. 8: (Color online) Phase diagram of the KLM of Eq. P, 
for Jhi ~ '2Jh2 = and n = 0.75. At weak coupling, 
< Jx/t ^ 1.3, the system is in a gapless dimer state with 
the small Fermi momentum kp = 7rn/2. A spin-gapped phase 
sets in at intermediate coupling 1.3 < J/c/t < 4, part of which 
retains nonzero dimer order (for 1.3 < JK/t < 1.7) and part of 
which where no dimer order is left (for 1.7 < Jx/i 4). The 
spin-gapped phase lacks a clear Fermi momentum. Finally, 
for 4 < Jx/t the system enters a paramagnetic LL phase 
with kp — Tv{n + l)/2 serving as the new Fermi momentum. 



FIG. 7: (Color online) Fourier transform of the density- 
density correlation function of Eq. ((Tj, for n — 0.75 and dif- 
ferent Kondo interactions. For Jk = there is a single sharp 
cusp at 2kp — 0.757r, in good agreement with the exact non- 
interacting result depicted by the red circles. For large Jk, 
the sharp structure at 2kp is replaced with two very shallow 
cusps at q — 2iT ~ 2kp — 0.257r and 2q, in agreement with 
C{k) for the corresponding t — J\ — J2 model^ (black dashed 
line). At the intermediate couplings JK/t ~ 2, 3, and 4 there 
is a sharp modulation at q, which degenerates into the shallow 
cusp of the t — J\ — J2 model upon increasing Jk- 



duce the exact noninteracting result for Jk — 0, and are 
consistent with the t — Ji — J2 model for large Jk- For 
Jk = there is a single sharp cusp at 2kp — 0.757r, which 
gradually smoothens as Jk is switched on. For large Jk 
(represented by J^/t — 16), there are two very shallow 
cusps at q and 2q, corresponding to modulations at 2kp 
and 4:kp- Both limits of small and large Jk are compat- 
ible with the Friedel oscillations observed in Fig. [51 

A more complex structure is found for the intermedi- 
ate couplings Jif/t = 2, 3 and 4. Here a sharp modula- 
tion develops at q, accompanied by additional wiggles at 
smaller values of k. We believe that the latter wiggles are 
a finite-size effect since they can be practically removed 
by using a windowing function that smoothly falls off to- 
ward the chain edges. By contrast, the primary peak at 
q not only persists but actually grows in magnitude if a 
windowing function is used. Similar to the Friedel oscil- 
lations, the primary peak at q evolves nonmonotonically 
with increasing Jk , first growing in magnitude before de- 
generating into the shallow cusp of the t — Ji — J2 model 
for large Jk- The overall behavior of C{k) in this regime 
is consistent with a tendency toward charge-density-wave 
or superconducting correlations that may accompany the 
opening of a spin gap. Lastly we note that C(fc) shows no 
discernable signature related to k* , reinforcing the spin 
origin of this modulation wavelength. 



IV. DISCUSSION AND CONCLUSIONS 

We begin by summarizing the zero-temperature phase 
diagram of the KLM of Eq. for J^ =2Jh2 = i/2 and 
n = 0.75. For weak Kondo couplings, < Jk/^ ^ 1-3, 
the system is in a gapless dimer state with the small 
Fermi momentum kp = 7rn/2 = 0.3757r. At intermediate 
Kondo couplings, 1.3 < Jk/^ ^ 4, a spin-gapped phase 
sets in,^^ part of which retains nonzero dimer order \D\ 
(for 1.3 < JK/t < 1.7), and part of which where no 
dimer order is left (for 1.7 < Jk/^ S 4). In the latter 
regime the system is presumably a Luther-Emery liquid. 
In contrast to the dimer state at weak coupling, the spin- 
gapped phase lacks a clear FS. There is no distinctive step 
in the electron momentum-distribution function, and the 
magnetic structure factor for the Kondo spins displays 
incommensurate modulations at 2k*, where k* is a new 
characteristic momentum. The latter wavelength is not 
associated with a new Fermi momentum, as neither the 
Friedel oscillations encoded in nr{k) nor the density- 
density correlation function C{k) show any discernable 
signatures related to fc*. Rather, above JK/t ~ 1.7 both 
nr{k) and C(fc) display pronounced modulations related 
to kp alone, in contrast to the regime JK/t 1-7 where 
nr{k) has simultaneous modulations at 2kp and 2k p- Fi- 
nally, for 4 < JK/t the system enters the paramagnetic 
LL phase of the effective Ji — J2 model, with kp serving 
as the new Fermi momentum. A sketch of the resulting 
phase diagram of the KLM is displayed in Fig. [51 

It should be emphasized that the opening of a spin gap 
in the KLM is by itself not new. A spin-gapped phase has 
long been establishec^ for the Kondo-Heisenberg model 
with n < 1 and Jh2 = 0, including recent indications 
for quasi-long-range superconducting correlations at a 
nonzero momentumi^ Here, however, the spin-gapped 
phase shows up as an intermediate phase, separating two 
paramagnetic phases with different Fermi momenta. This 
is quite different from previous results for Jh2 — Or 
where a single gapless phase was reported for values of 
Jk that exceeded the spin-gapped phase. 

Our results should be compared to those of Pivovarov 
and Si,^^ who used a perturbative renormalization-group 
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(RG) analysis to study a more general form of the KLM, 
including the effect of spin-exchange anisotropy in J^-^ 
and Jh'2 ■ For the SU(2) spin-symmetric interactions 
considered here, the perturbative RG and DMRG re- 
sults differ in two respects: (i) Perturbative RG pre- 
dicts a region of coexistence between the weak-coupling 
spin-Peierls phase and the strong-coupling Kondo-singlet 
phase, whereas no such region is found by the DMRG. (ii) 
The strong-coupling phase is predicted by RG to have a 
spin gap, while our DMRG study finds a LL. Indeed, we 
have verified by explicit calculations that the t — Ji — J2 
Hamiltonian onto which the system maps at strong cou- 
pling is gapless for rihoio = 0.75, thus forming a LL. Note, 
however, that a spin-gapped phase is expected in the 
i — Ji — J2 model as well when tuned sufficiently close to 
half filling (i.e., for n < I in the KLM). 

The transition studied in this work does not appear to 
be consistent with the local picture for the QPT in heavy 
fermion compounds fi^iii^ We did not observe a sudden 
change in the size of the FS, despite the fact that such 
a change would be more favorable in ID where the FS 
is reduced to just two isolated points. The intermediate 
spin-gapped phase for 1.3 < Jx/t < 4 seems to be a 
region where the Fermi surface reconstructs from kf. to 
kp. 

It is not obvious to what extent can the QPT studied 
in this work be compared to transitions involving true 
magnetic order, as seen experimentally in heavy fermions 
systems. For instance, the opening of a spin gap in the 
intermediate phase appears to be related to our choice 
of the Majumdar- Ghosh ground state for the isolated 
spin chain. Nevertheless, frustration can arise directly 
from the RKKY interaction itself, as is known to occur 
at quarter fiUingJ^ We have studied the Hamiltonian of 
Eq. ID) for n = 0.5 both with and without explicit frus- 
tration and found the same qualitative behavior. The 
electron momentum-distribution function n(k) displayed 
the same general behavior as for n = 0.75, including 
an intermediate region with no apparent Fermi momen- 
tum. At the same time, we did not find any spin gap 
in the intermediate regime, nor did we identify any new 
characteristic momentum k* associated with a shifting 



structure in either n{k) or S{k). The absence of k* for 
n = 0.5 likely stems from the fact that 2k f. and 2k p are 
indistinguishable for this particular filling factor. 

This naturally raises the question of how generic is 
the transition observed for n = 0.75. Is it representa- 
tive of other filling factors in Eq. ([1]) or does the QPT 
vary qualitatively as a function of 71? Preliminary results 
for n — 0.25 and n = 0.875 suggest the foUowingi^ For 
all filling factors studied there is an intermediate region 
where kp cannot be defined. This aspect, as well as the 
overall behavior of the momentum-distribution function, 
appears to be generic. However, details of the intermedi- 
ate region do depend on n. As stated above, we did not 
find any spin gap for n = 0.5, in contrast to the filling fac- 
tors n — 0.75 and 0.875 which are both spin gapped and 
in qualitative agreement with each other. The picture for 
n = 0.25 seems to be more complex, and may potentially 
involve more than one intermediate phase. This possi- 
bility, currently under study, may suggest a qualitative 
difference between the regimes n < 0.5 and 0.5 < n < 1. 

To conclude, we conducted an exhaustive investigation 
of the transition from a spin-Peierls phase with a small 
Fermi momentum to a LL phase with a large Fermi mo- 
mentum in a ID KLM. Our findings indicate a rather 
complex transition that exceeds the predictions of local 
criticality. It remains to be seen which of our results ex- 
tend to other electronic fillings and other variants of the 
ID KLM, let alone to higher dimensions. Further inves- 
tigations of this fascinating issue are clearly in order. 
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